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Abstract 

We perform a measurement of the branching fractions for t —> lyvv, {I - e,ii) decays for a minimum photon energy 
of 10 MeV in the r rest frame using 430 fb“^ of collisions collected at the center-of-mass energy of the T(45) 
resonance with the BaBar detector at the PEP-II storage rings. We find S(t fiyvv) - (3.69 + 0.03 + 0.10) x 10“^ and 
S(t ^ eyvv) = (1.847 + 0.015 + 0.052) x 10“^ where the first quoted error is statistical and the second is systematic. 
These results represent a substantial improvement with respect to existing measurements for both channels. 
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Leptonic t decays are generally well suited to inves¬ 
tigate the Lorentz structure of electroweak interactions 
in a model-independent way (ll] . Leptonic radiative de¬ 
cays T ^ lyvv, (/ = e, fj), in particular, have been stud¬ 
ied for a long time ||2l because they are sensitive to the 
anomalous magnetic moment of the t i|3i . Currently, the 
branching fraction of the t —> eyvv decay has been mea¬ 
sured only by the CLLO collaboration 101 using 4.68 
fb * of e^e^ collisions. Lor a minimum photon energy 
Ey,mm =10 MeV in the t rest frame they quote the re¬ 
sult (1.75 + 0.06 + 0.17) X 10“^, where the first error is 
statistical and the second systematic. The CLLO col¬ 
laboration also made the most precise branching frac¬ 
tion measurement of r —» fiyvv decay for a minimum 
photon energy in the t rest frame Ey^mm - 10 MeV, 
(3.61 + 0.16 + 0.35) X 10 ^ In addition, the OPAL col¬ 
laboration found S(t —> yuyvv) = (3.0 + 0.4 + 0.5)x 10“^ 
for a minimum photon energy of 20 MeV in the t rest 
frame 10 . In the present work, we perform a measure¬ 
ment of T —> lyvv, branching fractions for a minimum 
photon energy of 10 MeV in the r rest frame. 
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This analysis uses data recorded by the BaBar detec¬ 
tor at the PLP-II asymmetric-energy e^e^ storage rings 
operated at the SLAC National Accelerator Laboratory. 
The data sample consists of 431 fb * of e^e~ collisions 
recorded at a/s = 10.58 GeV/c. The expected cross sec¬ 
tion for T-pair production is cr^r = 0.919 + 0.003 nb 0 
corresponding to a data sample of about 400 million r- 
pairs. A detailed description of the BaBar detector is 
given elsewhere 101. Lor this analysis, a Monte Carlo 
(MC) simulation has been used to estimate the signal 
efficiency and to optimize the search. Simulated r-pair 
events are generated using KK2f 0 and t decays are 
simulated with Tauola llOtl . Linal-state radiative ef¬ 
fects in T decays are simulated using Photos d. A 
r-pair MC sample is generated where each r lepton de¬ 
cays to a mode based on current experimental knowl¬ 
edge d. A separate r- pair MC sample is generated 
where one of the r leptons decays to r ^ lyvv, and 
the other decays according to known decay modes. We 
exclude signal events in the former sample to obtain a 
T-pair background sample. In this article, charge con¬ 
jugation is always assumed. The MC simulated back¬ 
grounds samples include qq (mm, dd, ss, cc), 

and BB (B - B^, B^) events, where i^iT events are 
generated by KK2f 191], qq events are generated using 
the JETSET generator lfl3ll while BB events are simu¬ 
lated with EVTGEN |ll4]. The detector response is simu¬ 
lated with GEANT4 |^|. Two-photon and Bhabha back¬ 
grounds are estimated directly from the data. 
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The signature for t lyvv decays is a charged par¬ 
ticle (“track”), identified either as an e or a /r and a 
neutral deposit in the calorimeter (“cluster”) such that 
mass and energy of the lepton-photon pair are compat¬ 
ible with that of the parent t lepton. Events with two 
well-reconstructed tracks and zero total charge are se¬ 
lected, where no track pair is consistent with being a 
photon conversion in the detector material. The polar 
angle of each track in the laboratory frame is required 
to be within the calorimeter acceptance range to ensure 
a good particle identification. For every track, we re¬ 
quire pt > 0.3 GeV/c and a missing transverse momen¬ 
tum of pT,miss > 0.5 in the event. All neutral clusters 
are required to have a minimum energy of 50 MeV. We 
also reject events with neutral clusters with less than 110 
MeV if they are closer than 25 cm from a track, where 
the distance is measured on the inner wall of the EMC. 

Each event is divided into hemispheres (signal- and 
tag- hemisphere) in the center-of-mass (CM) frame by a 
plane perpendicular to the thrust axis, calculated using 
all reconstructed charged and neutral particles iQ.The 
signal hemisphere must contain exactly one track and 
one neutral cluster. The tag hemisphere must contain 
exactly one charged track, identified either as an elec¬ 
tron, muon or pion, and possibly one additional neutral 
cluster or (n =1, 2). Each is built up from a 
pair of neutral clusters with invariant mass consistent 
with that of the To suppress e^e~ and 

Bhabha events, we reject events in which the signal-side 
and tag-side leptons have the same flavor. In the signal 
hemisphere, we require for both channels that the dis¬ 
tance between the charged track and the neutral cluster 
on the inner wall of the EMC is less than 100 cm. For 
every event, the magnitude of the thrust is required to be 
within 0.9 and 0.995. The lower limit on the thrust mag¬ 
nitude helps to reject most qq events while the upper 
limit helps to reject e^e~ P^P and Bhabha events. 
For the same reason, we impose the total reconstructed 
energy to be less than 9 GeV. 

Electrons are identified applying a multivariate al¬ 
gorithm using as input the ratio of calorimeter en¬ 
ergy to the magnitude of the vector momentum of the 
track (E / p), the ionization loss in the tracking system 
(dEldx), and the shape of the shower in the calorime¬ 
ter. Muon identification makes use of a bagged decision 
tree (BDT) algorithm ifT^ . which uses as input the num¬ 
ber of hits in the IFR, the number of interaction lengths 
traversed, and the energy deposition in the calorimeter. 
Since muons with momenta less than 500 MeV/c do not 
penetrate into the IFR, the BDT uses also information 
obtained from the inner trackers to maintain a very low 
7T - p misidentification probability with high selection 


efficiencies. The electron and muon identification effi¬ 
ciencies are 91% and 62%, respectively. The probability 
for a ;7r to be misidentified as an e in r decays is below 
0.1%, while the probability to be misidentified as a p is 
around 1% depending on momentum. 

After the preselection, both samples are dominated 
by background events. For the r ^ pyvv sample, the 
main background sources are r ^ pvv, t —> jijf'v de¬ 
cays, —> p^p^ events, and t ^ ttv decays. For 

the T ^ eyvv sample, almost all background contri¬ 
bution is from r ^ evv decays in which the electron 
radiates a photon in the magnetic field of the detector 
(bremsstrahlung). Further background suppression is 
obtained by placing requirements on the angle between 
the lepton and photon in the CM frame (cos Oiy). For 
T — » pyvv we require cos Oiy > 0.99, while for t — > eyvv 
we require cos Oy > 0.97 (see Figs. [T]and|2]l. To reject 
background from t —> evv decays in the t —> eyvv sam¬ 
ple, we further impose a minimum value for the invari¬ 
ant mass of the lepton-photon pair My > 0.14 GeV/c^ 
for this channel. 


xl0^ 



Figure 1: Cosine of the angle between the lepton and photon momenta 
in the CM frame for radiative r decay into a muon after applying all 
selection criteria except the one on the plotted quantity. The selection 
criteria on the plotted quantity are highlighted by the vertical lines; we 
retain the regions indicated by the horizontal arrow. In green t —* /ivv 
decays, in blue r —> decays, in red —> lEyT events, in 

white signal r —» fiyvv decays, and in yellow other (background) r. 
The black dots are data. 

In addition to the afore-mentioned quantities, the se¬ 
lection criteria in both channels use the energy of the 
photon and the distance between the track and neutral 
deposit dy on the inner EMC wall. The selection cri¬ 
teria are optimized in order to give the smallest statis¬ 
tical and systematic uncertainty on the branching frac¬ 
tions. After optimization, for t —> pyvv, we require 
cos0/y > 0.99, 0.10 < Ey < 2.5 GeV, 6 < dy < 30 
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Figure 2: Invariant mass of the lepton photon pairfor radiative r decay 
into an electron after applying all selection criteria except the one on 
the plotted quantity. The selection criteria on the plotted quantity are 
highlighted by the vertical lines; we retain the regions indicated by the 
horizontal arrow. In green t —> ^vv decays, in white signal r —^ ^yvv 
decays, and in blue other (background) r. The black dots are data. 


Mode 

T ^ fuyyy 

T ^ eyyy 

Efficiency (%) 

0.480 + 0.010 

0.105 + 0.003 

BIN 

0.102 + 0.002 

0.156 + 0.003 

Nexp 

15649 + 125 

18115+135 

^obs 

15688 

18149 


Table 1: Signal efficiency, background contribution BjN = 
^bkgli^sig +Nbkg)^ where Nsig are signal and Nbkg background events, 
number of expected events {Nexp = Nsig+N^kg) and number observed 
events (A^obs) for the two decay mode after applying all selection cri¬ 
teria. 


cm and Miy < 0.25 GeV/c^. For the r ^ eyvv chan¬ 
nel, we require cos 0/y > 0.97, 0.22 < Ey < 2.0 GeV, 
^ < diy < 65 cm while the lower cut on the invariant 
mass is set to M/y > 0.14 GeV/c^. The selection ef- 
hciency determined using the MC samples is given in 
Table □ 

The branching fraction is determined using 


Elobs V fylig 

Icttt Ze 


( 1 ) 


where Nobs is the number of observed events, Nbkg is the 
number of expected background events, (Ttt is the cross 
section for r pair production, Z is the total integrated lu¬ 
minosity and e is the signal efficiency determined from 
the MC. After applying all selection criteria we hnd 


S(t ^ luyvv) = (3.69 + 0.03 + 0.10) x 10^^ (2) 


S(t eyvv) = (1.847 + 0.015 + 0.052) x 10^2(3) 


where the hrst error is statistical and the second is sys¬ 
tematic. Efficiency, background expectations (Nbkg) and 
the number of observed events (Vobs) are shown in Ta¬ 
ble [U 

Uncertainties on signal efficiency estimation and on 
the number of the expected background events affect the 
final result. For background estimation, we define con¬ 
trol regions that are enhanced with background events. 
For r —> fJ.yvv, where the major background contri¬ 
bution is not peaking in cos we invert the cut on 
cosOfjy. For cosOfjy < 0.8, the maximum expected sig¬ 
nal rate is 3% of the corresponding background rate. 
The maximum discrepancy between the MC sample 
prediction and the number of observed events is 8%, 
with an excess of events in the MC sample. We take 
this discrepancy as estimate of the uncertainty on back¬ 
ground prediction. For t —> eyvv whose major back¬ 
ground contributions have similar cos d^y distributions 
as signal, we apply a similar strategy after requiring the 
invariant mass M/y < 0.14 GeV/c^; in this case we take 
cos Oey < 0.90. The maximum contamination of sig¬ 
nal events in this region is 10%, and the maximum dis¬ 
crepancy between the prediction and the number of ob¬ 
served events is 4% with an excess of data events. We 
take this value as an estimate of the uncertainty on back¬ 
ground rate. The error on the branching fractions due 
to the uncertainty on background estimates are 0.9% 
for T —> iJ.yvv, and 0.7% for r —> eyvv, respectively 
(Table |2]l. Cross-checks of the background estimation 
are performed by considering the number of events ex¬ 
pected and observed in different sideband regions im¬ 
mediately neighboring the signal region for each decay 
mode and found to be compatible with the aforemen¬ 
tioned systematic uncertainties. 

The most important contributions to the error on ef- 
hciency come from the uncertainties on particle iden- 
tihcation and photon detection efficiency. Uncertain¬ 
ties on particle identihcation efficiency are estimated on 
data control samples, by measuring the variation of the 
data and MC efficiencies for tracks with the same kine¬ 
matic properties. The uncertainty on the efficiency of 
the electron identihcation is evaluated using a control 
sample consisting of radiative and non-radiative Bhabha 
events, while the uncertainty for muons is estimated us¬ 
ing an e^e^ —> i-i^n^y control sample. The uncertainty 
on the pion misidentihcation probability, as muon or 
electron, is evaluated using samples of t — » nnnv de¬ 
cays. The corresponding systematic error on the effi¬ 
ciency for T —> lyvv is 1.5% for both channels. To 
estimate the uncertainty on photon detection efficiency, 
we rely on two different processes depending on photon 
energy; for high energy photons we use e^e~ l^iTy 
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T juyvv 

T eyvv 

Selection Criteria 

- 

2.0 

Photon efficiency 

1.8 

1.8 

Particle Identification 

1.5 

1.5 

Background Evaluation 

0.9 

0.7 

PDG BF 

0.7 

0.7 

Nrr T pairs 

0.6 

0.6 

MC Statistics 

0.5 

0.6 

Trigger Selection 

0.5 

0.6 

Track Reconstruction 

0.3 

0.3 

Total; 

2.8 

3.4 


Table 2: Summary of systematic contributions to the branching frac¬ 
tion (in relative percent) for the two signal channels. 

events while for low energy photons we extract the un¬ 
certainty from 7r“ reconstruction efficiency, in which the 
photon kinematics can be fully reconstructed using the 
muon pair. Using fully reconstructed e^e~ —> 
events, data and MC are found to be compatible within 
1% for photon energies above 1 GeV. For photon en¬ 
ergies below 1 GeV we measure the reconstruction 
efficiency from the ratio of the branching fractions for 
T ^ nv and t ^ pv decays. The resulting uncertainty 
on the reconstruction efficiency is found to be below 
3%. Including the 1.1% uncertainty on the branching 
fractions, the resulting uncertainty on single photon de¬ 
tection efficiency is 1.8%. We take this last value as sys¬ 
tematic contribution to the efficiency for r —» lyvv due 
to photon detection efficiency. Another possible source 
of systematic uncertainty arises from the choice of the 
selection criteria; for t —> eyvv we observe a maximum 
deviation of 2% from the mean value of the branching 
fraction depending on the value of cos Oiy used for selec¬ 
tion. We take this value as estimate for the uncertainty 
on the efficiency due to the choice of selection crite¬ 
ria. A similar study on r —> pyvv shows that in this 
case there is no dependence of the result from selection 
criteria and thus the corresponding uncertainty is neg¬ 
ligible. All other sources of uncertainty in the signal 
efficiency are found to be smaller than 1.0%, including 
limited MC statistics, track momentum resolution, ob¬ 
servables used in the selection criteria, and knowledge 
of the tau branching fractions. 

In conclusion, we made a measurement of the branch¬ 
ing fractions of the radiative leptonic r decays t ^ eyvv 
and T —> pyvv for a minimum photon energy of 10 
MeV in the t rest frame using the full dataset of 
collisions collected by BaBar at the center-of-mass en¬ 
ergy of the T(45) resonance. We find S(t —> pyvv) - 
(3.69 + 0.03 + 0.10) X 10-^ and eyvv) = (1.847 + 


0.015 + 0.052) X 10“^ where the first error is statistical 
and the second is systematic. These results represent an 
improvement of about a factor of three for both chan¬ 
nels with respect to the previous experimental bounds 
0], reducing both statistic and systematic contributions 
for both channels. The main contribution to the to¬ 
tal error for both measurements comes from the pho¬ 
ton detection efficiency and particle identification. For 
T ^ eyvv, there is also an important contribution com¬ 
ing from the dependence of the final result from the se¬ 
lection criteria on the lower value of the outgoing elec¬ 
tron and the photon. Our results are in agreement with 
the SM values, S(t —> pyvv) - (3.686 + 0.009) x 10“^ 
and S(t —» eyvv) = 0.843 + 0.002) x 10“^, obtained 
from the TAUOLA HI MC, which uses the PHOTOS 
O package to simulate QED radiative corrections in 
the decay up to third order in the fine structure constant 
a. 
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